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■ Abstract. Recent determinations of the internal composition and structure 
^ , of two helium-atmosphere variable white dwarf stars, GD 358 and CBS 114, 

have led to conflicting implied rates for the ^■^C(a,7)^^0 reaction. If we 

7—i • assume that both stars were formed through single-star evolution, then the 

^ I initial analyses of their pulsation frequencies must have differed in some 

^ ■ systematic way. I present improved fits to the two sets of pulsation data, 

Q . helping to resolve the tension between the initial results. 
OO ■ 

o : 

O;! ■ 1. Introduction 

o ■ 

. When a white dwarf is being formed in the core of a red giant star during 

Qh| helium burning, there are two nuclear reactions that compete for the avail- 

Q ■ able helium nuclei: the 3a reaction, which combines three helium nuclei to 

^ , form carbon, and the ^^C(a,7)^^0 reaction, which combines an additional 

^ I helium nucleus with the carbon to form oxygen. At a given core tempera- 

■ ture and density, the relative rates of these two reactions largely determines 
I the C/0 ratio in the resulting white dwarf star. The rate of the 3a reaction 

^ ' is known to about 10% precision, but the ^^C(a,7)^^0 reaction is still un- 

certain by about 40%. So, if we can measure the C/0 ratio in the core of a 
pulsating white dwarf, it is effectively a measurement of the ^^C(a,7)^^0 
reaction rate. The C/0 ratio is interesting by itself, since the core composi- 
tion in our models affect the derived cooling ages of white dwarfs by up to 
a few Gyr Q. But we can also use it to provide an independent measure- 
ment of a nuclear reaction that is important to many areas of astrophysics, 
from the energetics of type la supernovae explosions to galactic chemical 
evolution. 

The model-fitting method that I describe below has only been applied to 
DB white dwarfs, since they are structurally the simplest. But in principle 
it can be extended to the pulsating DA stars quite easily, and with a little 
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more work to the DOVs. Presently, the method requires that the spherical 
degree of the pulsation modes is known, and sufficient data exist for only two 
stars — though we have just finished a Whole Earth Telescope run on a third 
object. The first application was to GD 358, which showed 11 consecutive 
radial overtones with the same spherical degree during a WET run in 1990 
m, |3|. The second application was only recently finished, and came from 
single-site data on the star CBS 114, which showed 7 independent modes 
that all appear to be £=1 Q. What we set out to do was search for a 
theoretical model that could reproduce, as closely as possible, the pulsation 
periods that we have observed in these stars. 

2. Model-Fitting 

We adjusted five different model parameters to try to match the observed 
periods. To make the final result as objective as possible, we wanted to 
explore the broadest range for each model parameter, defining the limits 
of the search based only on the physics of the model and on observational 
constraints. We allowed the mass to be anywhere from 0.45 to 0.95 M©, 
which encompasses the vast majority of known white dwarf masses |5|. The 
temperature of our models was allowed to vary from 20,000 to 30,000 K, 
which easily includes the spectroscopic temperature determinations of all of 
the known DBV white dwarfs whether or not trace amounts of hydrogen are 
included in their atmospheres Q . We looked at helium layer masses ranging 
from a fractional mass of lO"'^ where helium burning will begin at the base 
of the envelope, down to a few times 10^^, close to the limit where our 
models no longer pulsate 0]. The final two parameters describe a simple 
C/0 profile that has a constant ratio {Xq) out to some fractional mass 
point (g), where it then decreases linearly in mass to zero oxygen at the 
95% mass point. The important features from the standpoint of pulsations 
are the central C/0 ratio and the location and slope of the composition 
transition. 

To get reasonable resolution, we allowed 100 values for each param- 
eter inside these limits, so there are 10^*^ possible combinations of these 
five model parameters. Even if we had 1000 of today's fastest processors, it 
would still take more than a year to calculate all of these models, so we used 
a slightly more clever method employing a genetic algorithm to explore the 
many possibilities. The way genetic algorithms work is initially like a Monte 
Carlo method, where we just take a random sample of parameter combi- 
nations. After this initial random sampling, a genetic algorithm explores 
new regions of the search space based on a sort of survival of the fittest 
scheme. By passing simulated data through this process, we can quantify 
how long we need to let it run to find the correct set of parameters most 
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of the time, and by running the entire process several times with different 
random initiahzation, we can ensure a very high probabihty of finding the 
globahy optimal set of model parameters. In the end, the method requires 
a few million models to be calculated, concentrated mostly around the re- 
gions of the search space that yield better than average matches to the 
observations. So in addition to the optimal set of model parameters, we 
also end up with a fairly decent map of the search space, which gives us 
some sense of the uniqueness of the final answer. 

To learn what the optimal values of the mass and central C/0 ratio say 
about the ^^C(q,7)^^0 reaction, we need additional models that evolve 
a star from the main sequence through the red giant phase and into a 
white dwarf Q. For a given white dwarf mass, there are several things in 
the models that can be adjusted to change the central C/0 ratio, but the 
ingredient that affects it the most is the ^^C(a,7)^^0 reaction. So, all we 
have to do is adjust this rate in the models until we end up with the mass 
and central C/0 ratio that matches the fit from the genetic algorithm. 

3. Results 

In the initial application of this method to the data from the two stars, 
we found that CD 358 implied a reaction rate significantly higher than 
the extrapolations from laboratory measurements Q, while CBS 114 was 
right in line with the expectations [Q. This led us to speculate that there 
might be some systematic error affecting our analysis of the two stars in 
different ways. There was another slight worry in the initial results: the 
masses and temperatures for both stars differed significantly from those 
inferred from spectroscopy. We thought this may have resulted from the 
use of slightly different mixing length parameters than the spectroscopic 
studies, so we repeated the fits using ML2/a=1.25 to see if the discrepancy 
would disappear. In doing so, we also realized that there was a systematic 
difference in the way we were analyzing the two stars: for CD 358 the 
pulsation modes were consecutive radial overtones, so we were using both 
the periods and the spacings between the periods to judge which models 
provided the best match. For CBS 114 there was a gap in the sequence of 
observed modes, so we only used the periods themselves. When we repeated 
the fit for CD 358, we used only the periods to determine the best fit, just 
as we had for CBS 114. 

The results of the new fits are shown in Table 1. Notice that both objects 
yield a rate for the ^^C(a, 7)^^0 reaction that is consistent with laboratory 
extrapolations {Ssoo = 200 ± 80, 0), but CD 358 still seems to be a bit 
high compared to CBS 114, so this is only part of the answer. Also note 
that switching to the mixing length parameters used in the spectroscopic 
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TABLE 1. Optimal ML2/a=1.25 Models 



Object 


Tefl 


M/Mq 


l0g(MHe/M.) 


Xo 


9 


o-p(s) 


S300 


GD 358 


21,300 


0.695 


-2.95 


0.69 


0.49 


1.11 


215 ± 20 


CBS Hi 


20. ."•){){) 


0.71.") 


-6.77 


().")8 






1()() ± 20 



analysis did not resolve the differences between the implied masses and 
temperatures: the two methods of inferring these parameters do not agree. 
Curiously, the optimal model of GD 358 has a thick helium layer, but for 
CBS 114 it is thin. Finally, both models show the transition point from a 
constant C/0 ratio near the same fractional mass, and this location does 
not favor convective overshoot. Probably the decisive tests of all of these 
puzzles will come as we apply this method to additional DBV white dwarfs. 

4. Conclusions 

By measuring the interior composition of pulsating white dwarfs, we can 
get precise measurements of the important ■'^^C(a, 7)^^0 reaction rate, and 
as our models improve we can be more and more sure that they are not 
only precise, but also accurate. The model-fitting tool that we have used, 
involving a genetic algorithm, is a very powerful way to explore large ranges 
of interesting physical parameters and find the globally optimal model to 
match the observations. And since it evaluates so many models along the 
way, it produces some good maps of the search space in the regions where it 
is most interesting. True, this method is still computationally intensive, but 
there is no getting around that if we want the global solution, and Linux 
clusters are getting cheaper and cheaper. Finally, in the future we hope to 
be able to say something more about the detailed shape of the C/0 profile 
all the way from the center to the surface. 
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